function Hv = H(a,b)
%  H Возвращает вектор фокусировки для нашей антенной системы 
%  Входные фргументы - азимут и угол места

[kx, ky, kz] = sph2cart(b, a, 1);
k = [kx; ky; kz];

r{1} = [1/4, -1/4, 0]; % координаты антенн. координаты вектора r считаем в длинах волн
r{2} = [-1/4, -1/4, 0];
r{3} = [-1/4, 1/4, 0];
r{4} = [1/4, 1/4, 0];

NAnt = length(r) ;% число антенн

phi = nan(NAnt,1);
for i = 1:NAnt
phi(i) = 2*pi*r{i}*k;

end

Hv = exp(1i*phi);
end